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I. 


INTRODUCTION 


This  quarterly  technical  summary  report  covers  the/ work  performed 
during  the  period  January  through  March  1967.  Part  of /the  work  report¬ 
ed  herein  was  performed  under  Contract  AF33  (657 ) -15919,  which  covered 
SDL  operations  for  January  and  February  1967.  Work  previously  completed 
or  currently  in  progress  is  mentioned  only  as  it  relates  to  analyses 
completed  during  this  reporting  period. 

Analyses  completed,  for  which  results  have  been  reported,  are  dis¬ 
cussed  in  Section  II  under  descriptive  headings.  Section  III  contains 
a  discussion  of  the  support  and  service  tasks  performed  for  in-house 
projects  and  for  other  VELA-UNIFORM  participants.  Appendix  A  is  a  list¬ 
ing  for  those  organizations  receiving  SDL  data  services  during  this 
period . 

II.  WORK  COMPLETED 


A.  Detection  of  Surface  Waves  From  Small  Events  at  Teleseismic 
Distances 


One  of  the  important  measures  of  source  mechanism  is  the  ex¬ 
citation  of  surface  waves.  This  encompasses  such  estimates  as  total 
energy,  radiation  pattern,  frequency  spectrum,  and  radiation  of  Love 
wave  energy  relative  to  Rayleigh  wave  energy.  At  one  time  or  another 
practically  all  of  these  estimates  have  been  applied  in  the  study  of 
large  magnitude  events.  However,  because  the  surface  waves  from  small 
events  (M  <  5)  usually  arrive  at  teleseismic  distances  with  amplitudes 
at  or  below  the  noise  level,  special  techniques  must  be  employed  to 
extract  the  desired  information.  The  purpose  of  this  analysis  was  to 
use  a  matched  filter  approach  for  detecting  small  amplitude  surface 
waves,  estimating  their  total  energy  content  and  determining  their 
radiation  pattern. 


Basically  the  matched  filter  approach  amounts  simply  to  search¬ 
ing  a  record  x(t)  for  a  known  wave  form  y(t).  In  particular,  a  least- 
squares  approach  was  used  and  the  following  decision  criteria  were 
chosen  to  determine  whether  or  not  a  signal  was  presents 


1.  A  relative  maximum  in  the  envelope  of  C(t)  must  fall  in 
the  expected  time  window  based  on  travel-time  information. 


This  value  of  C  (t)  must  be  greater  than  or  at  least 
max 


comparable  to  typical  peak  values  in  the  envelope  of  C (t)  outside  the 
expected  time  window. 


Mmm 


3.  An  arrival  must  be  detected  within  the  proper  expected  time 
window  at  each  of  several  stations. 
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The  matched  filter  approach  requires  that  one  knows  the  de¬ 
sired  or  expected  waveform  in  order  to  search  effectively  for  that 
waveform  in  a  noisy  record.  In  order  to  obtain  a  suitable  "expected" 
waveform  and  at  the  same  time  assure  that  propagation  effects  would 
be  properly  accounted  for,  we  chose  the  surface  wave  from  a  larger 
event  in  the  source  region  of  interest  as  the  filter  y(t).  Thus,  no 
matter  how  complicated  the  paths  of  propagation  w^re,  the  paths  for 
other  events  in  the  same  source  region  would  be  nearly  coincident  with 
those  of  the  larger  event  which  produced  y(t),  so  the  only  major  dif¬ 
ferences  in  recorded  waveform  would  be  source  differences „ 

In  order  to  test  the  effectiveness  and  sensitivity  of  the 
matched  filter  approach  for  surface  waves,  several  test  cases  were 
investigated,  being  adding  a  known  signal  to  actual  noise  at  different 
aignal-to-rtoise  ratios,  effects  of  varying  epicentral  distance, 

and  the  effects  of  amplitude  spectra  on  matched  filter  response. 

The  matched  filter  approach  was  then  applied  to  actual  data. 
Results  for  three  different  source  regions  -  Aleutians  (Amchitka) , 
Kamchatka,  and  Hawaii  -  are  given  below.  The  observation  stations  are 
LKSM  stations  or  VELA  observatories  in  North  America.  The  epicenter 
information  for  .the  three  sets  of  events  are  given  in  Table  I. 

The  organization  in  the  figures  to  follow  conforms  to  the 
standard  format  adopted  for  displaying  all  matched  filter  results 
for  actual  events.  That  is,  a  station  set  consists  of  four  traces? 
the  top  trace  is  the  signal  y(t),  the  second  trace  is  the  observed 
seismogram  x(t),  the  third  trace  is  C  (t),  and  the  bottom  trace  is 
C (t ) .  The  legend  at  the  left  gives  values  of  distance,  a,  and  CC 
as  well  as  station  and  event  identification.  The  arrows  indicate 
the  location  of  the  expected  signal  arrival  at  each  station.  This 
is  where  the  peak  in  the  matched  filter  output  should  occur?  that  is, 
the  beginning  point  of  y(t)  in  x(t).  The  interval  between  each  time 
mark  is  50  seconds  unless  otherwise  indicated. 

1.  Aleutian  Source  Region 

The  LONGSHOT  seismograms  were  scanned  for  surface  waves  with  the 
earthquake  signal  used  as  the  matched  filter.  Results  showed  that 
the -Rayleigh  waves  were  clearly  detected  at  all  but  one  station. 

Figure  1  is  an  example  of  these  results.  The  Rayleigh  wave  amplitude 
of  LQNGSHQT  relative  to  that  of  the  magnitude  5.9  earthquake,  however, 
was  below  .12  everywhere  in  the  distance  range  analyzed  (3500-7000  km) . 
This  means  that  the  surface  wave  magnitude  for  LONGSHOT  was  only  about 
5.0  if  the  surface  wave  magnitude  of  the  reference  earthquake  was  equi¬ 
valent  to  its  body  wave  magnitude.  It  also  means  that  the  R/P  exci¬ 
tation  for  the  explosion  was  nearly  an  order  of  magnitude  smaller  than 


Table  I  -  USCGS  Study  Epicenter  Information 
Used  for  Matched  Filter 


Figure  1.  Matched  Filter  Analysis  o 
Rayleigh  Waves  with  22  November  1965 
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the  R/P  excitation  for  an  equivalent  size  earthquake  - 

Figure  2  shows  the  matched  filter  results  for  the  magnitude  4.0 
earthquake  (20  February  1966)  listed  in  Table  I.  Notice  that  the 
source  depth  for  this  event  was  48  km  which  means  that  the  surface 
wave  excitation  is  necessarily  smaller  than  if  the  same  source  were 
located  closer  to  the  surface.  From  Figure  2  it  appears  that  this 
event  was  barely  detected  at  several  stations  but  not  detected  at 
every  station.  The  values  of  4  average  about  .01  which  is  approx¬ 
imately  the  value  expected  for  a  magnitude  difference  of  2.  Certain¬ 
ly  this  event  is  at  or  near  the  lower  limit  of  detectability  using 
the  matched  filter  approach. 

2 .  Kamchatka  Source  Region 


From  Table  I,  it  can  be  seen  that  in  these  cases  the  observing 
distances  were  in  the  range  5,300-8,000  km  which  overlaps,  but  ex¬ 
tends  further  than  the  distances  investigated  for  the  Aleutian  events. 
It  can  also  be  seen  that  all  these  events  are  located  within  2-3 
degrees  of  each  other. 

Results  from  the  5.3,  5.1,  and  5.0  magnitude  events  clearly 
showed  that  the  events  were  positively  detected.  Shown  in  Figure  3 
are  some  of  the  matched  filter  results  for  the  21  January  1965  4.4 
magnitude  event .  It  should  be  noted  that  this  was  a  deep  focus  event 
(119  km) ,  hence  equivalent  to  a  smaller  magnitude  shallow  event  as 
far  as  Rayleigh  wave  excitation  is  concerned.  It  appears  that  this 
event  was  barely  detected  by  the  method  with  estimates  of  a  around  .03. 
This  value  gives  a  magnitude  estimate  of  about  4.2  for  the  21  January 
event,  which  seems  large,  considering  the  source  depth  of  119  km. 

In  any  case  this  event  appears  to  define  the  threshold  of  detectability 
for  the  Kamchatka  region. 


Hawaiian  Source  Region 


-  Results  from  both  the  4.4  and  4.1  magnitude  events  shows  that 
both  were  apparently  detected.  In  figure  4  are  shown  the  results  for 
the  shallow  focus  (11  km)  magnitude  4.1  event  of  13  August  1964.  This 
event  is  apparently  detected  at  most  of  the  stations.  However,  the 
corresponding  correlation  coefficients  and  the  estimates  of  a  are  small. 
The  apparent  magnitude  of  the  small  event  is  about  1.6  smaller  than  the 
reference  event.  However,  the  filtering  effect  of  oceanic  propagation 
of  Rayleigh  waves  will  greatly  attenuate  periods  shorter  than  about 
15  seconds  and  as  the  peak  in  excitation  shifts  to  shorter  periods  for 
the  smaller  magnitude  events,  this  effect  should  cause  very  signifi¬ 
cant  reductions  in  surface  wave  energy,  hence  reduce  the  matched  filter. 
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Figure  2.  Search  for  Rayleigh  Waves  From  a  Magnitude 
4.0  Event  at  Depth  40  km.  Using  22  November  1965 
Magnitude  5.9  Event. 
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Figure  3.  Search  for  Rayleigh  Waves  From  a  Magnitude  4.4 
Event  at  Depth  119  km.  Using  the  Magnitude  5.8  Event  of  29 

January. 1965. 
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It  has  been  demonstrated  experimentally  that  the  surface  waves 
from  small  events  can  be  detected  at  teleseismic  distances  using  a 
matched  filter  method  and  that  reasonable  estimates  of  some  source 
characteristics  can  be  made.  What  has  not  been  demonstrated  is  the 
minimum  threshold  of  detection  for  the  method.,  The  underlying  dif¬ 
ficulty  which  always  can  affect  the  estimated  threshold  seriously  is 
that  the  body  wave  magnitudes  used  as  an  obsolute  standard  for  com¬ 
paring  events  may  be  badly  in  error .  A  more  reliable  absolute  measure 
of  source  strength  is  needed  to  evaluate  the  performance  of  the  matched 
filter  in  detecting  actual  events „ 

One  of  the  most  useful  features  of  the  matched  filter  method  is 
that  it  permits  one  to  make  meaningful  surface  wave  magnitude  esti¬ 
mates  even  when  the  signal  to  noise  level  ,  is  so  low  that  one  cannot 
visually  detect  the  presence  of  the  signal.  This  means  that  the 
matched  filter  method  can  be  used  to  obtain  reliable  surface  wave  magni¬ 
tude  estimates  for  events  so  small  that  meaningful  AR  estimates  of 
magnitude  cannot  be  obtained. 

It  is  concluded  that. 


1.  In  principle  the  matched  filter  approach  can  be  useds 

(a)  to  detect  weak  surface  wave  signals  for  signal  to  noise 
ratios  as  low  as  about  0.35; 

(b)  to  make  reliable  relative  magnitude  estimates  for  S/N  as 
low  as  .5; 

(c)  to  determine  radiation  pattern  relative  to  a, reference 
event ; 

(d)  to  estimate  the  general  shape  of  a  small  event's  am¬ 
plitude  spectrum  relative  to  that  of  the  reference  event. 

2.  The  results  presented  for  a  number  of  actual  events  have  dem¬ 
onstrated  the  practical  applicability  and  usefulness  of  the  matched 
filter  approach  in  studying  weak  teleseismic  surface  waves. 

3.  The  matched  filter  used  in  conjunction  with  the  surface  wave 
arrival  time  patterns  appropriate  for  different  source  regions  of 
interest  allows  one  to  discriminate  against  all  but  those  events  in 
the  particular  source  region  of  interest. 

4.  Rayleigh  waves  produced  by  the  LONGSHOT  explosion  were, de¬ 
tected  at  teleseismic  distances.  The  excitation  of  LONGSHOT  Rayleigh 
waves,  however,  was  an  order  of  magnitude  smaller  than  the  Rayleigh 

wave  excitation  for  an  equivalent  magnitude  earthquake  in  the  same 
region „ 

. 

5.  Love  waves  could  not  be  detected  for  the 'LONGSHOT , explosion . 


Teleseismic  surface  wave  signals  could.be  detected  for  small 


itifs 


different  source  regions  investigated.  The  actual  lower  limit  for 
each  region  is  still  uncertain. 


The  (pP-P)  Time  Difference 

mere  is  a  difference  in  the  travel-time  anomaly  for  P  and  pp 
at  the  TFO  Extended  Array.  This  difference  is  sufficient  in  some 
cases  to  result  in  near  cancellation  of  pP  if  the  seismograms  are 
summed  by  aligning  on  P.  On  most  seismograms  it  is  not  possible  to 
determine  the  time  of  pP  with  sufficient  accuracy  to  justify  the 
above  statement.  If,  however,  the  examples  are  limited  to  seis¬ 
mograms  with  well-defined  pP,  then  the.  scatter  of  the  time  difference 
(pP-P)  can  be  demonstrated. 

Four  groups  of  seismograms  were  selected  from  the  vertical 
seismograms  recorded  at  the  TFO  Extended  Array  for  the  following 
earthquakes  % 

Chile-Bolivia  Earthquake  -  16  April  1965 
Peru  Earthquake  -  10  May  1965 
J.  Coast  Chile  Earthquake  -  20  July  1965 
Kurile  Is.  Earthquake  -  18  May  1965 

„Tbe  *lme  difference  (pP-P)  observed  at  the  9  stations  of  the 
TFO  Extended. Array  for  the  above  earthquakes  were  plotted  on  Figure 
'  gainst  A  (surface  distance).  In  each  case  the  scatter  of  the 

def?nLW«S  ^bout  °”5  seconds.  The  points  were  read  by  timing  a  we_- 
defined  peak  or  through,  rather  than  attempting  to  read  first  motion. 

~ m J!he •  PmrJ?et lines  in  Figure  5  were  plotted  from  the  Jeffreys- 
Shxmshoni  Tables  (Jeffreys  and  Shimshoni  1964).  They  establish  a  Y 

me^uSinSi.WhiCh  the  SlOPe  and  SGatter  of  the  pP-P  times  may  be 

0  01  a^d'o  Parame!riC  yalUa°  of  VR  for  values  other  than  0.00, 

0.01  ad  0.02  were  interpolated  from  the  tables. 

3  well-defined  phases  and  selected  earthquakes 

the  times  for  both  P  and  pP  are  believed  to  be  correct  with  O.l  second. 

feren«e(npfpras'aUil0nt-able^  pradict  an  “=tease  in  the  time  dif- 

measured  bv  h/n  w^  £unftlon  °f  *  (surface  distance)  and  depth  of  focus 

eartVs  radius  18  hypocentral  aeS*h  minus  33  km.  and  R  is  the 

earth  s  radius.  This  increase  is  measured  by  the  slope  of  the  h/H  Dints 

of  Figure  5.  it  varies  from  0.1  second  for  the  KurUe  Earthquake  to 
tt9  „1™sthr  °-fSeCOndS  for  the  Chile-Bolivia i Earthquake^acfoss 

swell  If  5®  dlstances  measured  by  the  TFO. Extended  Array.  The 

pread  of  surface  distances,  observed  by  the  TFO  Array,  is  commolly 

3  to  4  ,  varying  with  azimuth. 
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'5h&  visual  evidence  from  the  scatter  of  the  plotted  v 
(pP-P)  is  clearly  shown.  In  order  to  obtain  a  statistical 
of  this  scatter  the  following  analysis  was  made.  The  rela 
smallness  of  the  sample  limits  the  conclusions  to  an  indie 
the  confidence  limits. 


The  mean  and  variances  of  the  6 
computed  usirg  the  standard  formulas 


The  95%  confidence  intervals  of  the 
from  the  following  formula 


If  the  observed  times  for  (pP-P)  are  all  0.1  second  greater  or 
less  than  the  mean,  the  time  differences  of  a  pair  becomes  0.2  seconds 
The  phase  difference  of  1.0  sec  period  pulse  is  1/5  x  360°  or  72°. 

This  would  be  sufficient  phase  difference  to  produce  some  distortions 
however,  the  signal  would  probably  be  recognizable.  If  this  assumptioi 
is  used  as  a  basis  for  computation  then 


0.0367  >  a  >  0.00457 

This  is  the  95%  confidence  region,  for  the  assumption  each  time  is 
0;1  sec  greater  or  less  than  the  mean  is  the  above  range;  This 
difference  (72°  for  1  sec  period)  can  be  thought  of  as  a  threshold. 

If  the  time  differences  are  greater,  say  0.2;  then  the  phase  shift 
is  144°  and  summing  produces  near  cancellation .  This  criteria  of 
cancellation  is  not  readily  est abl ished ,  s ince  the  pulses  to  be 
summed  contain  a  range  of  frequencies/  It  is  probable  that  f or  re¬ 
inforcement  to  occur,  when  the  pulses  are  summed,  requires  a  narrower 
range  of  variance  than  the  range  computed  for  the  0.1  second  assumpt- 
i.0n«  In  each  cnee  the  actual  95%  confidence  interval  is  greater  than 
that  computed  for  the  0.1  sec.  assumption* 
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The  two  filters  selected  for  evaluation  of  performance  were: 

1.  2-Pole  Phaseless  High  Pass  Filter  -  this  is  set  with  the 
3  db  point  at  1  cps  and  only  suppresses  the  .2  cps  peak  in  the 
noise.  It  gave  no  visible  distortion  of  any  of  the  signals,  including 
first  motion. 

2.  Second  Order  Gapped  Finite  Difference  operator  with  the  first 
spectral  peak  at  1.1  cps  and  the  first  null  at  2.2  cps.  It  suppresses 
the  .2  cps  and  the  2  cps  noise  peak  and  gave  only  slight  distortion 

in  the  1  cps  signal  but  gave  considerable  distortion  in  the  1.8  cps 
signal. 


Noise  was  obtained  from  a  pseudo  random  number  generating  function, 
tested  and  shown  statistically  to  provide  a  sequence  of  uncorrelated 
numbers.  A  synthetic  signal  was  obtained  by  taking  the  displacement 
potential  as  an  impulse  into  a  system  which  models  linear  creep  by 
a  dashpot  and  spring,  in  series. 

This  noise  and  signal  were  used  to  evaluate  these  two  filters. 
Results  confirmed  earlier  findings  from  using  real  data.  The  gapped 
secord  difference  filter  appears  to  yield  the  best  performance  as  a 
simple  detector.  The  indicated  gain  in  S/N  ratio  of  8.5  db  and  gain 
in  threshold  of  7.5  db  are  compared  with  the  high  pass  filter's 
performance  of  3.5  db  and  2.5  db,  respectively.  One  of  the  keys  to 
optimum  selection  of  bandpass  filters  appears  to  involve  the  sup¬ 
pression  of  spectral  peaks  occurring  in  the  noise  spectrum,  especial¬ 
ly  the  0.2  cps  and  2.0  cps  peaks.  How  this  is  best  done  with  minimum 
signal  distortion  is  not  resolved  here:  however,  the  gapped  difference 
operator  was  the  best  design  evaluated.  Any  distortion  introduced  by 
the  filter  can  be  effectively  removed  after  array  processing. 

^ •  Noise  Analysis  of  Single  Channel  Deqhostinq  Filters 


The  purpose  of  applying  a  single  channel  inverse  deghosting  operator 
is  to  suppress  the  echo  reflected  at  the  earth's  surface.  ~For  teleseismic 
signals,  the  incidence  angle  is  approximately  15°,  thus  the  down-going 
pulse  cou.ld  be  modeled  with  a  reflection  coefficient  of  0.9  and  the 
up-going  pulse  with  a  reflection  coefficient  of  1.1.  In  ^practice,  however, 
reflection  coefficients  greater  than  1.0  lead  to  unstable  operators,  so 
that  single  channel  echo-suppression  is  feasinle  oniy  if  it  is  attempted 
to  partially  remove  the  echo  in  isolating  the  up-going  and  down-goinq 
pulses. 
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Two  different  deghosting  methods,  two  different  bandpass  filters 
and  unfiltered  data  were  the  five  cases  evaluated  by  performance  arrays. 
To  produce  each  element  in  the  array  we  operated  on  four  minutes  of 
data °  each  array,  16  hours  of  data,  and  for  all  five  cases,  80  hours 
of  data.  The  fact  that  all  of  these  performances  figures  were  computed 
With  less  than  40  minutes  of  computer  time  attests  to  speed  of  the  band¬ 
pass  and  deghosting  operators. 

Results  showed  that  arrays  of  performance  values  showing  .the  change 
m  noise  variance  due  to  deghosting  indicate  significant  sensitivity 
to  both  the  deghosting  method  and  pre-band  pass  filtering,  ..As  a  spe¬ 
cification  for  single  channel  deghosting,  the  apparent  reflection  coef¬ 
ficient  used  to  derive  the  operator,  should  be  constrained  such  that 
the  noise  amplification  does  not  exceed  a  specified  design  criteria 
such  as  1  db„  The  apparent  reflection  coefficient  used  to  design  the 
echo  suppressor  is  observed  as  highly  dependent  on  the  echo  time, 

III .  SUPPORT  AND  SERVICE  TASKS 


A,  VELA-UNIFORM  Data  Services 

As  part  of  the  contract  work-statement,  the  SDL  provided 
one  or  more  of  the  following  support  and  service  functions  for  VSC 
and  other' VELA  participants s 

-  copies  of  16  and  35  mm  film 

-  playouts  of  earthquakes  and  special  events 

-  copies  of  existing  composite  analog  tapes 

-  composite  analog  tapes  of  special  events 

-  use  of  1604  computer  for  checking  out  new, programs 

-  or  running  production  programs 

-  copies  of  digital  programs 

-  digitized  data-in  standard  formats  or  special  formats 
for  use  on  computers  other  than  the  1604 

-  running  SDL  production  programs,  such  as  power  spectral 
density  and  avray  processing  on  specified- data 

-  digital  x-y  plots  of  power  spectra  or  digitized  data 

-  signal  reproduction. booklets 

-  copies  of  MIT  Geophysics  Program  Set  II 

-  space  for  visiting  scientists  utilizing, SDL  facilities 

to  study  data  and  exchange  information  with  SDL  personnel 


During  this  report  period,  52  such  projects  were  com¬ 
pleted  and  the  14  organizations  receiving  these  services  are  listed 
in  Appendix  A,  ^  ,v  v 


B.  Data  Library 


The  Data  Library  contains  digitized  seismograms,  digital  com¬ 
puter  programs  and  composite  analog  magnetic  tapes,  all  available 
for  use  by  the  VELA-UNIFORM  program. 

The  following  additions  were  made  during  this  period. 

1.  Digital  Seismograms  -  317  including 

-  data  from  32  explosions  and  5  underwater  events 

-  noise  samples  from  LASA,  TFO,  UBSO ,  CPSO,  and  WMSO 

-  noise  from  several  vertical  arrays 

-  earthquake  data  recorded  at  various  stations  from  15  events 

2.  LASA  Data  -  127  digital  tapes 

-  in  addition  to  the  LASA  digital  tapes  in  the  library,  there 
is  also  a  master  calibration  tape  which  contains  the  magni¬ 
fication  (digital  counts  per  millimicron)  of  each  sensor 
for  every  subarray.  These  magnifications  have  been  computed 
for  all  calibration  tapes  currently  in  house.  As. each  new 
calibration  is  received,  it  is  routinely  run  through  the 
program  CALIBR  and  added  to  the  master  tape. 

3.  Digital  Programs  -  JLl  including: 

APCTTAT  - 

To  take  an  input  tape  containing  a  label  and  (NC)  no.  of 
data  channels,  obtain  all  possible  combinations  of  NC  taken  two 
at  a  time,  and  form  an  output  tape  containing  a  label  with  each 
set  of  two  records. 

.  COHERENCY  - 

To  compute  all  the  ordinary  coherencies,  auto-spectra,  and 
phase  relations  of  a  set  of  input  data  channels. 

FKSPTRUM  - 

To  compute  and  display  the  frequency-wave  number  power  spectra 
of  seismic  noise  along  with  a  response  function  for  the  correspond¬ 
ing  array.  This  program  is  intended  .as  a  replacement  for  certain 
parts  of  program,  PEAKAY  using  the  Cooley-Tukey  method  to  estimate 
auto  and  cross  power  spectral  density  functions. 

FSTMOFIL  - 

To  perform  recurrsive  filtering  of  seiamic  data  using  two 
methods,  namely  RECFIL  and  SDL,  in  order  to  determine  the  one 
which  preserves  first  motion  best. 


MAKESUC  - 

To  read  an  input  seismo 
containing  the  altered  labe 
may  be  detrended,  scaled,  s< 
label  is  changed  accordingl; 


MAXLIKE-21  - 

This  Fortran-63  program 
maximum- likelihood  digital  : 
to  seismic  array  data. 


RODRDATA  - 

To  take  an  input  tape  cc 
records  (a  label  and  so  man] 
generate  a  new  tape  with  an] 
any  set  re-written  such  as  1 


SIMVARRY . - 

To  simulate  an  N  company 
different  criteria. 


TWX  - 

To  compute  a  least-mean- 
one  channel  of  an  array  fron 
apply  the  filter  to  construe 
the  error  time  function,  anc 
error. 


VBRPROC 


XMNG 


GREELEY 

STERLING 

CHASE  V 

PRESTERLING 

One  event  for  Geotec 

Three  events  for  Tex 
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APPENDIX  A 


ORGANIZATIONS  RECEIVING  SDL  DATA  SERVICES 
JANUARY  -  MARCH  1967 


Dominion  Observatory 
Earth  Sciences,  Teledyne 
General  Atronics  Corporation 
Geotech,  Teledyne 
Oregon  State  University 
Pennsylvania-  State  University 
Princeton  University 
RADC 

Southern  Methodist 1  University 
'St.  Louis . University 
Texas  Instruments,  Inc. 
Underwater  Systems,  Inc. 

U...  S.  Coast  and  Geodetic  Survey 
Vitro  Corporation 
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